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We determine the ratio F K jF n in QCD with Nf — 2 + 1 flavors of sea quarks, based on a series 
of lattice calculations with three different couplings, large volumes and a simulated pion mass 
reaching down to about 190 Me V. We obtain Fk/F k = 1.192(7) aaf (6) SVJ with all the sources of 
systematic uncertainty under control. 
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1. Introduction 

We determine Fk/F k in QCD through a series of dynamical lattice calculations such that all the 
sources of systematic uncertainty are properly taken into account [jlj|. We use Nf = 2 + 1 dynamical 
quarks, with a single quark whose mass is close to the physical strange quark mass (m s ~ mf ys ), 
and two degenerate flavours of light quarks heavier than in the real world u and d quarks, but 
with masses varying trough a range that allows a controlled extrapolation to the physical point. 
Concerning finite volume effects, the spatial size L is large enough so that the values of F^/F^in 
our ensembles can be corrected for small finite-volume effects. We simulate at three different 
values of /3 to have full control over the continuum extrapolation. 

2. Simulation and analysis details 

Here we will not give any details about our actions for the gauge and fermion fields, or about 
the algorithms used for the simulation. The interested reader should consult [||]. 

To set the scale and adjust the quark masses, we use aM n , aM^ and either aM? or aM^ (to 
estimate the systematics, as we will see). We extrapolate for each value of the lattice spacing the 
values aM n ,aMK,aMz to the point where any two of the ratios agree with the experimental values. 
We subtract electromagnetic and isosping breaking effects to the experimental values of the masses: 
we use M£ hys = 135MeV, Atf ys = 495 MeV and M| hys = 1318MeV with an error of a few MeV, 
according to [Q]. 

Regarding m lu i, we have pion masses in the range 190 — 460 MeV. Fk/F k is measured with 
the valence quark masses equal to the sea quark masses (no partial quenching). 

The same dataset was successfully used to determine the light hadron spectrum 

3. Treatment of the theoretical errors 

3.1 Extrapolation to the physical mass point 

We simulate with a strange quark mass already close to its physical value, but this is not the 
case for the light quarks. Thus our values of Fk/F k need an extrapolation to the physical point. 
There are three possible guides for this extrapolation: SU (3) chiral perturbation theory, heavy kaon 
SU (2) chiral perturbation theory, and Taylor fits. 

For the case of SU (3) chiral perturbation theory we use the expression for the ratio at NLO 
as a function of the measured pion and kaon masses. It is important to note here that the apparent 
convergence of Chiral perturbation theory is a statement that depends both on the observable and 
the statistical accuracy of the data. With our data, the ratio Fk / F n is well described by the NLO 
expression, but this does not mean that NLO expressions can be used in general to describe others 
quantities (see [§]). The ratio Fk/F„ is probably a benevolent quantity as chiral logs in F n and Fx 
cancel in part. 

SU (2) chiral perturbation theory [||] and its heavy kaon variant [{J determines the functional 
dependence of the ratio of decay constants on pion mass. We find that the ratio of NLO expressions 
from these two references describes well the data in our ensembles. 
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Having data in the range 190 MeV — 460 MeV, it is natural to consider an expansion about a 
regular point which encompasses both the lattice results and the physical point [Q, gj]. The mass 
dependence of Fk/F k in our ensemble and at the physical point is well described by a low order 
polynomial in M 2 and M\. 

Since all the three frameworks describe well our data, we will use all of them in our analysis, 
and use the difference between them to estimate the systematic uncertainty. 

3.2 Continuum limit 

Fk/Fji is an SU (3)-flavor breaking ratio, so that cutoff effects must be proportional to m s — 
m u d, that guided by SU (3) chiral perturbation theory we may substitute with M\ — M\. 

Cutoff effects are both theoretically (they partially cancel in the ratio) and in practice small 
numerically. In |§] we found that although our action is only formally improved up to order a, 
these small cutoff effects seem to scale like a 2 up to about 0.16 fm. Even if this is the case we 
can not exclude the possibility of cutoff effects proportional to a. Thus we have considered the 
following three options compatible with our data: no cutoff effects, and a flavor breaking term 
proportional to a ox a 2 . 



3.3 Infinite volume limit 

Stable states in a box with periodic boundary conditions have different masses and decay 
constants than the corresponding states in infinite volume. The difference vanishes exponentially 
fast with the mass of the lightest state in the box [JSj] . In our case, masses and decay constants are 
corrected with terms proportional to exp(— M n L). In our simulations M n L^A making finite volume 
corrections small. Moreover the sign of leading correction in F n and Fk is the same, so that they 
partially cancel in the ratio. 



Within chiral perturbation theory, the 1-loop [Q, [TOp and 2-loop [JTTp corrections for the pion 
and kaon decay constants are known, so we have decided to correct the values of our simulations 
with the 2-loop expression before fitting the data. The 1-loop expression will be used to estimate 
the systematic uncertainty due to finite volume corrections (see below). 



Similarly, we also correct the meson masses [11]. 



4. Fitting strategy and treatment of theoretical errors 

Our goal is to obtain Fk/F^ at the physical point, in the continuum and in infinite volume. To 
this end we perform a global fit which simultaneously extrapolates or interpolates M 2 — > M^| p h ys , 
M\ — > M\ I phys and a — > 0, after the data have been corrected for very small finite volume effects 
using the two-loop chiral perturbation theory results discussed above. To assess the various sys- 
tematic uncertainties associated with our analysis, we perform a large number of alternative fits. 

Excited states contribute to the correlators, so to estimate their effect, we choose 18 different 
time intervals dominated by the ground state. The difference between masses and decay constants 
coming from different time intervals are used to estimate the uncertainty associated with excited 
states. 

Scale setting systematic uncertainty is estimated by using M n , Mk and either Ms or Mq. 
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The chiral extrapolation systematic error is estimated in two ways. First we consider two 
different ranges of pion masses for the fits: 190 — 350 MeV and 190 — 460 MeV. Second, we 
consider a total of 7 different functional forms to extrapolate to the physical point. 3 of them 
come from the NLO SU (3) chiral perturbation expression: the ratio of decay constants (cancelling 
terms proportional to La), a NLO expansion of the ratio and a similar NLO expansion but for the 
inverse ratio F k /Fk- 2 functional forms come from the SU (2) chiral perturbation theory expression: 
the expanded ratio of decay constants, and a similar expression but for F k /Fk- The last 2 forms 
correspond to a Taylor fit for Fk/F k , and a Taylor fit for F k /Fk- It is important to note that the 
difference between these 7 functional forms provides an estimate of both higher order contributions 
within a concrete framework (e.g NNLO terms in SU (3) chiral perturbation theory), and systematic 
bias coming from a particular framework. 

As discussed above, cutoff effects are parametrised in three different ways: we consider fits 
with and without 0(a 2 ) and 0(a) corrections, as described in Sec. 3.2 . 

Following this procedure we have 18x2x2x7x3 = 1512 global fits. One of the 1512 
fits (corresponding to a specific choice for the time intervals used in fitting the correlators, scale 
setting, pion mass range, . . . ) can be seen in (Fig. [Taj). We emphasize that the % 2 per d.o.f. for our 
correlated fits are close to one. 

Although all these methods seem to describe well our data, not all of them do it in the same 
way, so we weight with the fit quality the central value of each fit. These 1512 weighted values 
can be used to construct a distribution, whose median is an estimate of the typical result of our 
analysis, therefore our desired final result (see fig. [T^). The width of the distribution is a measure 
of the systematic error of our analysis, thus we take the 16-th/84-th percentiles as an estimate of 
the systematic error of our computation. 

Finite volume effects are treated separately because we know a priori that the two-loop ex- 
pressions of [|ll]] are the most accurate expressions available and they describe well these effects in 
our data. To estimate the error associated with the finite volume effects, we repeat the full analysis 
using the 1 loop expression to correct the ratio Fk / F n and we also repeat the full analysis correct- 
ing only the value of F % (this can be seen as an upper bound to the real correction in Fk/F k ). The 
weighted (with the quality of the fit) standard deviation of these three values is used as an estima- 
tion of the uncertainty due to finite volume effect, and added to our systematic error by quadratures 
to produce the final systematic error. 

To determine the statistical error, the whole procedure is bootstrapped with 2000 samples, and 
the standard deviations of the 2000 medians used as our estimate of the statistical error. Our final 
error is computed as the sum by quadratures of the total systematic and statistical errors. 



5. Results 



Our final result for the ratio of decay constants is 



Fk 



1.192(7) stat (6) syst or 



Fn 



phys 



0.839(5) stat (4) syst 



(5.1) 



phys 



at the physical point, where all sources of systematic error have been included. 
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(a) Chiral extrapolation of the lattice data to the physi- 
cal point for a particular choice of time interval and mass 
cut (M n < 460 MeV in this case). Here we use a Taylor 
ansatze with no cutoff effects. The m s dependence has 
been subtracted to plot the data as a function only of M K . 




(b) Distribution of values of F^jF K . The large background 
distribution (yellow) represents the values of Fj^/F n ob- 
tained with different extrapolation formulas, pion mass 
cuts, parameterization of cutoff effects, time intervals and 
different methods to set the scale. Also shown is the fi- 
nal error interval (dashed lines) and the final value (black 
solid vertical line). 



Figure 1: And example of mass extrapolation and the distribution of fits used to obtain the final 
result and systematic error. 



Figure |a| shows our final result compared with the determination of F^/F^ from other dynam- 
ical lattice computations. There are two Nf = 2 computations by JLQCD [|l3|] and ETM [14]. With 
2+1 fermion flavours, we have a number of results obtained using MILC configurations: MILC 
[§, 0], NPLQCD @, HPQCD/UKQCD |Q and Aubin et al. iff]. The results by RBC/UKQCD 
[f7|] and PACS-CS [ 19] were also obtained with Nf = 2 + 1 simulations but with different configu- 
rations. It is worth noting that these results show a good overall consistency when one excludes the 
outlier point of [13]. 



6. Contributions to the systematic error 

Having estimated the total systematic error, it is interesting to decompose it into its individual 
contributions. To quantify these contributions, we construct a distribution for each of the possible 
alternative procedures corresponding to the source of theoretical error under investigation. These 
distributions are constructed by varying over all of the other procedures and weighing the results 
by the total fit quality. Then, we take the weighted standard deviations of the medians of these 
distributions as an estimate of the systematic uncertainty associated with the source of error under 
consideration. 



Table (|2b|) shows the estimation of the different sources of systematic error in our computation. 
Even having pion masses down to 190 MeV, the chiral extrapolation remains the main source of 
systematic error. This error is broken in two parts by changing the fit range, and using different 
expressions to extrapolate the data. In the Fig. || we can see how fits corresponding to different 
pion mass cuts, and different functional forms contribute to the final distribution of values. The 
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(a) Comparison of recent lattice computations of (b) Breakdown of the total systematic error on F^/F^ 
Fk/F%. into its various components, in order of decreasing im- 

portance. 

Figure 2: In the figure, we can see a comparison between our result and recent unquenched com- 
putations of Fjc/Fjz. The table shows, in order of importance, the different sources of systematic 
error. 





(a) Final distribution of values and the contribution to this (b) Final distribution of values and the contribution to this 
corresponding to different pion mass cuts. corresponding to different functional forms. 



Figure 3: Analysis of chiral extrapolation error. 



medians of the small constituent distributions are used to compute the error associated with each 
source of error, as was mentioned before. 

The next source of systematic error in importance are cutoff effects. As was commented 
earlier, we have three possible parametrizations for the cutoff effects: no cutoff effects at all, a 
flavor breaking term proportional to a, and a flavor breaking term proportional to a 2 . In Fig. pjajthe 
corresponding distributions are shown. In the same way Fig. |5a| shows the contributions coming 
from the two possibilities for setting the scale. 

The two remaining sources of theoretical error deserve a separate comment. First, the con- 
tamination with excited states is studied by using a total of 18 different fitting ranges for the cor- 
relators, corresponding to t m \ n /a = 5 or 6, for j8 = 3.3; 7, 8 or 9, for j8 = 3.57; 10, 11 or 12 for 
j6 = 3.7. In Fig. ^ we can see the final distribution compared with the distributions corresponding 
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(a) Final distribution of values and the contribution to this (b) Final distribution of values and the contribution to this 
corresponding to different cutoff effects. corresponding to two different fitting ranges (of the 18 

possibilities). The total area of the distributions have been 
rescaled to make them visible. 

Figure 4: Distributions with its contributions coming from different sources of theoretical error. 
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(a) Final distribution of values and the contribution to this (b) Final distribution of values, compared with the dis- 
corresponding to different scale settings. tributions corresponding to a different (1-loop and upper 

bound) finite volume correction. 



Figure 5 : Distributions with its contributions coming from different sources of theoretical error. 



to t m [ n /a = 5,7, 10 and t^ m /a = 6,9, 12. These distributions have been rescaled (so that they add 
to the total area of our final distribution), to make the small distributions more visible. 

Second, only the 2-loop finite volume corrections are included in the 1512 fits used to deter- 
mine the final central value, because we do not want to bias this value with the 1-loop contributions, 
that are, a priori, less accurate than the 2-loop corrections. To estimate the uncertainty associated 
with the subtraction of finite volume effects, we repeated the full analysis with 1-loop finite vol- 
ume corrections and with an upper bound to the correction, computed by including only the 2-loop 
correction to F % . Our final distribution of values with these two alternatives are plotted in Fig. |5^. 
The error associated with the finite volume effects is computed as the weighted (by fit quality) 
standard deviation of the medians of these three distributions, and added by quadratures to the 68% 
confidence interval of our final distribution to produce the final systematic error. 

One final comment about the procedure to obtain the final systematic error. The addition of 
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different procedures to construct the final distribution only can increase the systematic error. For 
example in figure < j3b| ) we can clearly see that dropping any of the procedures to extrapolate to the 
physical point (for example not using SU (3) fits), gives a final value well in our final error band, but 
a smaller systematic error. This general statement remains true for the other sources of systematic 
error. 
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